Low field phase diagram of spin-Hall effect in the mesoscopic regime 
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When a mesoscopic two dimensional four-terminal Hall cross-bar with Rashba and/or Dresselhaus 
spin-orbit interaction (SOI) is subjected to a perpendicular uniform magnetic field B, both integer 
quantum Hall effect (IQHE) and mesoscopic spin-Hall effect (MSHE) may exist when disorder 
strength W in the sample is weak. We have calculated the low field "phase diagram" of MSHE in 
the (B, W) plane for disordered samples in the IQHE regime. For weak disorder, MSHE conductance 
GsH and its fluctuations rms(GsH) vanish identically on even numbered IQHE plateaus, they have 
finite values on those odd numbered plateaus induced by SOI, and they have values Gsh = 1/2 and 
rms(GsH) = on those odd numbered plateaus induced by Zeeman energy. For moderate disorder, 
the system crosses over into a regime where both G s h and rms(GsH) are finite. A larger disorder 
drives the system into a chaotic regime where G s h = while rms(GsH) is finite. Finally at large 
disorder both G s h and rms(GsH) vanish. We present the physics behind this "phase diagram". 

PACS numbers: 71.70.Ej, 72.15.Rn, 72.25.-b 



Many recent papers have been devoted to the physics 
of spin-Hall effect 1] and a particular focus is the intrin- 
sic spin-Hall generated in non-magnetic samples by spin- 
orbital interaction (SOI) [2|, yj • So far, several experimen- 
tal papers have reported observations of spin-Hall effect 
in compound semiconductors and other systems 4]. The- 
oretically, it has been shown that for two dimensional 
(2D) samples in the clean limit, the Rashba SOI gener- 
ates a spin-Hall conductivity having a universal value of 
e/87r 3]. The presence of weak disorder destroys spin-Hall 
effect in large samples 0,0- In particular, a consensus ap- 
pears to have been reached in the literature that spin-Hall 
effect in disordered samples generated by linear Rashba 
SOI vanishes at the thermodynamical limit 0,0,0- 

For mesoscopic samples, numerical studies have pro- 
vided evidence that the mesoscopic spin-Hall effect 
(MSHE) can survive weak disorder |1 [H E3- For a 
four-probe disordered sample, MSHE conductance Gsh 
and its fluctuations rms(GsH) have been calculated for 
both linear Rashba and Dresselhaus SO interactions [ToL 
H^ . It was found [l^ that when the system is in the diffu- 
sive regime, the fluctuations rms(GsH) take a universal 
value with the same order of magnitude as the average 
Gsh itself, and is independent of the system size L, the 
disorder strength W the electron Fermi energy and the 
SO interaction strength. 

The situation becomes very interesting and more com- 
plicated when a perpendicular uniform external magnetic 
field B is applied to the 2D sample [l^. In this case, 
Gsh and rms(GsH) become functions of B. Most im- 
portantly, a magnetic field B can produce edge-states 
which are responsible for the integer quantum Hall effect 
(IQHE). Similar to the well known studies of the global 
phase diagram of quantum Hall effect [3, it will be very 
useful to map out the low field "phase diagram" of MSHE 
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FIG. 1: Transmission coefficient T12 versus E or B(T). For 
setup-I: (a) and (b). For setup-II: (c) and (d). Inset of (a): 
schematic plot of the setup-I; inset of (b): the corresponding 
flow of edge states. Inset of (c): schematic plot of the setup-II; 
inset of (d): the corresponding flow of edge states. 



in terms of the field strength B and the disorder strength 
W. Such a diagram allows one to clearly understand the 
role played by the edge-states and disorder. It is the 
purpose of this work to present this MSHE "phase dia- 
gram" for four-probe 2D disordered mesoscopic samples 
with linear Rashba and/or Dresselhaus SO interactions. 

Here we put "phase diagram" in quotes because the 
physics we study is mesoscopic, namely for samples in 
the coherent diffusive regime characterized by the rela- 
tion between relevant length scales, I < L < £. Here L 
is the linear sample size, I the elastic mean free path and 
£ the phase coherence length. As such, the "phases" 
in the "phase diagram" are states with zero or finite 
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values of Gsh and rms(GsH), an d no phase transi- 
tions are implied between these states. In particular, 
we found that with low disorder when IQHE is well es- 
tablished, both Gsh and rms(GsH) are zero identically 
on the even numbered IQHE plateaus, while they take 
finite values on the SOI dominant odd numbered IQHE 
plateaus. For Zeeman dominant odd numbered IQHE 
plateaus, G s h = 1/2 and its fluctuation vanishes. As the 
disorder is increased, both Gsh and rms(GsH) become 
nonzero when any edge-state is destroyed by the disorder 
in any IQHE plateau. Further increase of disorder brings 
the system to a "chaotic" regime where Gsh — while 
rms(GsH) ^ 0, finally at even larger disorder both Gsh 
and rms(GsH) vanish. These behaviors are organized in 
the low field phase diagram which we determine in the 
rest of the paper. 

We consider a 2D four-probe device schematically 
shown in the inset of FigHt (call it setup-II). A MSHE 
conductance Gsh is measured |lQ] across probes labeled 
2,4 when a small voltage bias is applied across probes 1 
to 3 so that a current flows between them. Gsh can be 
measured the same way when there is a uniform exter- 
nal magnetic field B which exists everywhere including 
inside the leads. Gsh is theoretically calculated from 
spin-current defined as I s = ft/ 2(1^ — I\) where I^i are 
contributions from the two spin channels. Note that the 
definition of I s is, in f act, in debate for regions where 
SO interaction exists |7l llq. To avoid this ambiguity we 
assume that in our device the SO interaction only exists 
in the shaded region (setup-II in FigH), namely in leads 
1,3 and in the central scattering region, but does not 
exist in leads 2,4 where we measure spin-current. This 
way, I s is well defined as above. For discussion purposes, 
we have also considered a device (setup-I, inset of FigGlO 
where SO interaction is present everywhere including in- 
side leads 2, 4. 

In the presence of linear Rashba interaction a so z • (a x 
k) with k = k + (e/hc)A, the Hamiltonian of the four- 
probe device is: 

H = y ] ^nm c nm a C nma + 9s ^ ] C nmcr( (7 ' ~B)aa' c nma' 
nma nmaa' 

~ t ^ ] i C n+l,ma C nma e ^ + C n,m-la C nma + h.C.] 
nma 

— t so [c n) r m j r i a {io~ x ) C7C7 >c nrnc7 > 

nmaa' 

~ cUl,ma( i( Ty)™' C nma>e- irn11 + h.C.} (0.1) 

where c^ ma is the creation operator for an electron with 
spin o~ on site (n,m), e nmcr = At is the on-site energy, 
t = h 2 /2/ia 2 is the hopping energy and t so = a so /2a is 
the effective Rashba spin-orbit coupling, g s = (l/2)g/iB 
(with g = 4) is the Lande g factor. Here n = Tiuo c /2t and 
uo c = eB/fic is the cyclotron frequency. Throughout this 
paper, we use t as the unit of energy. For L = 40a = 1/im, 
t = 1.5 x 10 _3 eV, and t so = 0.2t corresponds to a so = 



9 x 10- 12 eV.m|3. We choose A = (-By, 0,0) so that 
the system has translational symmetry along x-direction 
(from lead 1 to lead 3). Static Anderson- type disorder 
is added to q with a uniform distribution in the interval 
[—W/2, W/2] where W characterizes the strength of the 
disorder. The spin Hall conductance G s h is calculated 
from the Landauer-Buttiker formula 9] 

G sH = (e/87r)[(T 2Til - T 2U1 ) - (T 2T>3 - T 2i>3 )] (0.2) 

where transmission coefficient is given by T2 a ,i = 
Tt^G^^). Here G r ' a are the retarded and ad- 
vanced Green's functions of central disordered region of 
the device which we evaluate numerically. The quantities 
are the line width functions describing coupling of the 
leads to the scattering region and are obtained by calcu- 
lating self-energies due to the semi-infinite leads using a 
transfer matrices method The spin-Hall conductance 

fluctuation is defined as rms(G s #) = \J {G 2 sH ) — (G s h) 2 \ 
where (• • •) denotes averaging over an ensemble of sam- 
ples with different disorder configurations of the same 
strength W. The devices in Fig^ have L x L central 
square, and without losing generality we fixed L = 40 
grid points in our numerics. 

Before presenting the numerically determined "phase 
diagram" for the physics of MSHE using setup-II, let's 
first discuss the general physics of spin-Hall current. For 
this purpose we use setup-I where the SOI is everywhere 
so that the discussion is simpler. We first examine the 
spin-Hall "phase diagram" in the absence of SOI. In a 
magnetic field, edge-states are formed. FigEt,b shows 
transmission coefficient T12 for setup-I, which measures 
the number of edge-states, versus Fermi energy E or mag- 
netic field B. We observe that T12, or the number of 
edge-states, increases as E for a fixed B and it decreases 
as B is increased for a fixed E. Notice that the number 
of edge-states N can be either even or odd. The odd TV 
region in E or B is very narrow and is due to the Zee- 
man splitting that breaks the spin degeneracy. When N 
is even, spin-Hall current vanishes because all the edge- 
states are fully polarized with half of them pointing to 
one direction (say spin-up) and the other half pointing 
to opposite direction (spin-down). When N is odd, the 
spin-Hall conductance is 1/2. At weak disorder when 
all the edge-states survive, we therefore conclude that 
G s h — when N is even and G s h = 1/2 when N is odd. 
Furthermore, it is useful to examine fluctuations of the 
spin-Hall conductance rms(G S H) for these edge-states: 
we expect no fluctuations for all edge-states. As disor- 
der strength W is increased, we reach a point where at 
least one of the edge-states is destroyed and the system 
is in a spin-Hall liquid state characterized by G s h 7^ 
and rms(G S H) 7^ for any N. Further increasing W, 
we expect strong scattering to bring the system into a 
chaotic state of MSHE, characterized by G s h = and 
rms(G S H) 7^ 0. At even larger W, the system enters a 
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spin-Hall insulator state where G s h = rms(G S H) = 0. 

Next, we turn on the SOI and discuss its effect on 
the "phase diagram" . Fig^,d show transmission coeffi- 
cient T\2 for setup-I versus E or B for a fixed Rashba 
SOI t so = 0.2. We observe that the behavior of T12 
is similar to that of FigEt,b except that the region of 
odd TV is now much larger. When N is even, spin-Hall 
current vanishes as before. In the region of B when N 
is odd, two cases occur due to the competition between 
SOI which tends to randomize the spin polarization and 
the Zeeman energy which favors spin polarization along 
a fixed direction. If Zeeman energy is large enough, then 
G s h = 1/2 as before with rms(G S H) — while if SOI 
dominates then there is at least one edge-state that has 
both spin- up and down components: our numerical re- 
sults show that the composition depends on systems pa- 
rameters. As a result, there is a net spin-Hall current 
when N is odd. This discussion becomes clearer when 
we examine setup-II where the spin direction can be de- 
fined. At weak disorder when all the edge-states survive, 
we have the same conclusion as before, i.e. G s h = 
when N is even and G s h 7^ when N is odd. We expect 
no fluctuations for even N and for those odd N edge- 
states with G s h = 1/2, but finite fluctuations for the 
rest of odd N edge-states. Hence, at weak disorder, we 
have a "phase" of edge-state induced spin-Hall insulator 
with even N characterized by G s h = rms(G S H) = 0; a 
"phase" of edge-state induced spin-Hall liquid (but fluc- 
tuationless and Zeeman dominant) with odd N charac- 
terized by G s h = 1/2 and rms(G S H) — 0; and finally 
a "phase" of edge-state induced spin-Hall liquid (SOI 
dominant) with odd N characterized by G s h 7^ and 
rms(G S H) 7^ 0. As we increase the disorder strength, 
the "phase diagram" evolves through three regimes simi- 
lar to the case when SOI is off: a spin-Hall liquid regime, 
a chaotic regime, and a spin-Hall insulating regime. 

The discussion in the last paragraph gives the entire 
expectation for the low field MSHE "phase diagram". 
The problem of this discussion is that the spin-Hall cur- 
rent is not well defined in regions where SO interaction 
exists 05 [HI such as setup-I of FigGJi. Therefore, in the 
rest of the work we consider setup-II where SO interac- 
tion does not exist in leads 2, 4 so that spin-Hall current 
is well defined and measurable without ambiguity. The 
extra complication of setup-II is that there is an inter- 
face between spatial region with t so = and that with 
t so 0. This interface acts as a potential barrier caus- 
ing additional scattering of edge-states. In particular, at 
certain energies one of the edge-states goes directly from 
lead 1 to lead 3 due to this interface scattering. Insets of 
Fig. lb and Fig. Id show schematically the edge states for 
setup-I and II, respectively. In the inset of Fig. Id, how- 
ever, an edge-state is now transmitted directly from lead 
1 to lead 3 due to the interface scattering just discussed. 
We have confirmed that this is a generic feature which 
occurs at different Fermi energies. For a fixed Fermi en- 
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FIG. 2: (color online) (a). The edge state plateaus in (B,W) 
plane, (b). The transmission coefficient T12 for setup-I, setup- 
II, as well as direct transmission coefficient T13 as a function of 
B in the absence of disorder, (c). The spin-Hall conductance 
in (B,W) plane, (d). The spin-Hall conductance fluctuation 
in (B, W) plane. 



ergy, this can also happen when B is varied. In Fig. 2b, 
we plot the T12 for setup-I, and T12, T13 for setup-II, at 
W = 0. We observe that N = odd edge-states are much 
easier to be scattered while the TV = even edge-states are 
stable against interface scattering. Therefore, the regions 
in the MSHE "phase diagram" where TV = even becomes 
larger for setup-II than for setup-I. For instance, the mag- 
netic field B for the onset of N = 2 edge-state changes 
from 1.32T to 1.2T due to the interface scattering (for 
a device with lead width L = l/im). We emphasis that 
except for this extra complication of interface scatter- 
ing in setup-II, the general physics discussion of MSHE 
"phase diagram" for setup-I in the last paragraph, holds 
perfectly for setup-II. 

Fig. 2a depicts numerical result for the number of edge- 
states N as we vary B and W. We observe that the 
edge-states are gradually destroyed from the subband 
edge (measured in the lead 1) to the subband center 
when W is increased. From Fig. 2a we also observe that 
N = 2 edge-states are more stable against disorder than 
that of N = 3. Fig. 2c, d show spin-Hall conductance 
and spin-Hall conductance fluctuation, respectively, for 
W < 4 [13. They are perfectly consistent with the gen- 
eral discussion given above, namely G s h and rms(G S H) 
are finite for N =odd edge-states and in regions when at 
least one edge-state is destroyed by disorder. 

Fig. 3 plots the main result of this work, the low field 
"phase diagram" of MSHE. In the numerical calculations 
of this "phase diagram", we have computed 61 values of 
B, 40 values of W from W = to W = 4, and for each 
pair of (B, W) we averaged over 1000 impurity config- 
urations. The integers in the "phase diagram" indicate 
the number of edge-states N. At weak disorder, there 
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FIG. 3: (color online) The low field "phase diagram" of meso- 
scopic spin-Hall effect in (B, W) plane. Note that for disorder 
strength between W — 2 to W — 20, the system is the spin- 
Hall liquid. 

are three possible states: the TV = even edge-state in- 
duced spin-Hall insulator, the SOI dominant TV = odd 
edge- state induced spin-Hall liquid state, and the Zeeman 
dominant N = odd edge-state induced fluctuationless 
spin-Hall liquid. Since large magnetic field favors Zee- 
man term, so in N = odd plateau the SOI dominant spin- 
Hall liquid appears first for low magnetic field and crosses 
over to Zeeman dominant fluctuationless spin-Hall liquid 
at higher field. As W increases, the edge-states become 
destroyed and the system enters spin-Hall liquid where 
G s h 7^ and rms(G S H) ^ 0. A chaotic state of MSHE 
with G s h = and rms(G S H) 7^ is reached when W is 
increased further. Finally, the system enters a spin-Hall 
insulator state where G s h = = rms(G S H) = at large 
enough disorder. While this "phase diagram" is obtained 
for a particular value of Rashba SO interaction t so , we 
have checked that the general topology is the same for 
other values. In addition, the MSHE "phase diagram" in 
the (t SOl W) plane for a fixed magnetic field has similar 
features. We have also determined the phase boundary 
between the chaotic state of MSHE and spin-Hall insula- 
tor that are shown in Fig. 3 with the same resolution[l^. 

We have so far focused on linear Rashba SOI. A sim- 
ilar analysis can be carried out for Dresselhaus SOI by 
adding a term /3 so (a x k x — (J y k y ) in Ea. (|().l[) . It is well 
known that in the absence of Zeeman energy one has 
Is H ( a so = 0,(3 SO ) = Is H (a S o,Pso = 0) and I z sH (a so = 
ft so) — 0. Therefore, in the absence of Zeeman energy, 
the MSHE "phase diagram" for Dresselhaus SOI is the 
same as that of the Rashba SOI. In the presence of Zee- 
man energy, our numerical results for Dresselhaus SOI 
give a similar "phase diagram" . When both Rashba and 
Dresselhaus terms are present, a similar "phase diagram" 



is also obtained numerically for t so = 0.2 and t SQ 2 = 0.4 
(tso2 = Pso/2a). 

In summary, we have determined the low field "phase 
diagram" of mesoscopic spin-Hall effect. The "phase di- 
agram" is characterized by values of G s h and rms(G S H) 
in the W) plane and the main features include a spin- 
Hall liquid behavior where both G s h and rms(G S H) are 
nonzero, and by spin-Hall insulator behavior where both 
quantities vanish. Furthermore, the spin-Hall liquid can 
be induced by N =odd edge- states in weak disorder, and 
by destroying edge-states for larger disorder. The spin- 
Hall insulator behavior, on the other hand, is induced by 
N =even edge- states, and by very large disorder. The 
MSHE "phase diagram" is found to be true for both lin- 
ear Rashba and Dresselhaus SO interactions. 
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